Luminosity distance in Swiss cheese cosmology with randomized voids and galaxy 

halos 



Eanna E. Flanagan, 1, 2 [j Naresh Kumar, ^[j and Ira Wasserman 1, 2 [^| 

1 Laboratory for Elementary Particle Physics, Cornell University, Ithaca, NY 14853, USA 

2 Center for Radiophysics and Space Research, Cornell University, Ithaca, NY 14853, USA 

We study the fluctuations in luminosity distance due to gravitational lensing produced both by 
galaxy halos and large scale voids. Voids are represented via a "Swiss cheese" model consisting 
of a ACDM Friedman-Robertson-Walker background in which a number of randomly distributed, 
spherical regions of comoving radius 35 Mpc are removed. A fraction of the removed mass is 
then placed on the shells of the spheres, in the form of randomly located halos, modeled with 
Navarro-Frenk-White profiles. The remaining mass is placed in the interior of the spheres, either 
smoothly distributed, or as randomly located halos. We compute the distribution of magnitude 
shifts using a variant of the method of Holz & Wald (1998), which includes the effect of lensing 
shear. In the two models we consider, the standard deviation of this distribution is 0.065 and 0.072 
magnitudes and the mean is -0.0010 and -0.0013 magnitudes, for voids of radius 35 Mpc, sources 
at redshift 1.5, with the voids chosen so that 90% of the mass is on the shell today. The standard 
deviation due to voids and halos is a factor ~ 3 larger than that due to 35 Mpc voids alone with 
a 1 Mpc shell thickness which we studied in our previous work. We also study the effect of the 
existence of evacuated voids, by comparing to a model where all the halos are randomly distributed 
in the interior of the sphere with none on its surface. This does not significantly change the variance 
but does significantly change the demagnification tail. To a good approximation, the variance of 
the distribution depends only on the mean column depth and concentration of halos and on the 
fraction of the mass density that is in the form of halos (as opposed to smoothly distributed): it is 
independent of how the halos are distributed in space. We derive an approximate analytic formula 
for the variance that agrees with our numerical results to < 20% out to z ~ 1.5. 



I. INTRODUCTION 
A. Background and Overview 

A number of surveys are being planned to determine 
luminosity distances to various different astronomical 
sources, and to use them to constrain properties of the 
dark energy or modification to gravity that drives the cos- 
mic acceleration. Perturbations to luminosity distances 
due to gravitational lensing by large scale and galaxy 
scale structures are a source of error for these studies, 
see, e.g., Refs [THS]- 

In this paper we use the computational method de- 
veloped in Ref. [T] to study the effect of density in- 
homogeneities on luminosity distances in two idealized 
"Swiss cheese" models (3H7HHI] of large scale (~ 30 Mpc) 
and galaxy scale structures. Our models seek to capture 
the property that most of the matter is concentrated in 
galaxy halos on the outer edges of voids while the void in- 
teriors are relatively sparse. Our first model is an exten- 
sion of our previous work [l\ where we idealize the interior 
of a spherical void as a uniform underdense region and 
the surface of the sphere contains a randomized distribu- 
tion of galaxy halos with Navarro-Frenk-White (NFW) 
profiles [T3] • Our second model retains the randomly 
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distributed galaxy halos on the surface of the voids, and 
replaces the interior uniform density with randomly dis- 
tributed galaxy halos with NFW profiles. In both models 
we keep fixed the parameters of the voids and halos. Even 
though neither of these models represent realistic matter 
distributions within a void, they should capture the main 
qualitative features of lensing. 



B. Summary of Computation and Results 

In Section II, we give a detailed description of the NFW 
halo density profile that we adopt which is motivated by 
observations. We also describe our two models of the 
entire distribution of matter, including both large scale 
void structures and smaller scale halo structures. 

In Section III, we review the method we use to com- 
pute the distribution of lensing magnifications. We then 
describe how to compute the lensing convergence for a 
single halo. We compare our numerical results for the 
distribution of magnifications with analytic expressions. 
We also estimate the number of realizations required to 
get a reasonable accuracy in the computed distribution of 
magnifications (e.g, obtain the mean of the distribution 
to an accuracy of < 1% ). The accuracy of the numer- 
ical results scales as N~ x / 2 as expected, where N is the 
number of realizations. 

In Section IV, we study our first Swiss cheese model. 
We start by describing the model, and study the propa- 
gation of light rays through just a single void. Then, we 
derive analytic results for the magnification and use these 
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to check our numerical results. We study the expected 
number of halo intersections, the redshift dependence of 
the magnification distribution, and determine the contri- 
bution of shear to our results. We note that the standard 
deviation is ~ 3 times larger than that due to voids with 
no halos, specifically the model of [1] consisting of voids 
of radius 35 Mpc with smooth underdense interiors and 
a smooth overdensity concentrated on the surface of the 
sphere with a thickness of 1 Mpc. We show that the 
redshift dependence of the mean and standard deviation 
agrees with analytic results to < 10%. We also note that 
the standard deviation changes by less than 3% if shear 
is neglected (see Section IV C below). 

One effect which our models do not include is the clus- 
tering of halos, that is, the correlations between the loca- 
tions of different halos. While it would be more realistic 
to include the effects of clustering, our simplified models 
should capture the essence of the effects of large scale 
inhomogeneities. 

In Section V, we study the second Swiss cheese model. 
Again, we first describe the model and study just a sin- 
gle void. We derive analytic results for the lensing con- 
vergences and use these to check our numerical results. 
There is a higher probability of demagnification; this shift 
is expected because the density contrast inside the void 
is now sharper because it is empty (with a smattering of 
a small number of halos) whereas the first model has a 
smooth interior matter distribution. The redshift depen- 
dence of the mean and standard deviation of the second 
model arc similar to those of the first model. In the two 
models we consider, the standard deviations of this dis- 
tribution are 0.065 and 0.072 magnitudes and the means 
are -0.0010 and -0.0013 magnitudes, for voids of radius 35 
Mpc, sources at redshift 1.5, with the voids chosen so that 
90% of the mass is on the shell today. We compare the 
distributions for configurations with and without voids 
for a source at z s = 1.5. We find that the voids do 
not significantly change the variance but do significantly 
change the demagnification tail and the mode. 

We find that since the distribution is skewed, the mode 
is positive, while the variance is determined primarily by 
rays that intersect halos. The scale of the voids does not 
significantly influence our results. The main parameters 
that determine the mode and variance of the distribu- 
tion is the mean column depth and concentration of ha- 
los and the fraction of the mass density that is in the 
form of halos (as opposed to smoothly distributed). The 
distribution of halos in space (i.e., in the interior versus 
the surface) is unimportant. Hence, our models bracket 
the range of possibilities of magnifications. Our analysis 
is generally consistent with other analytic and compu- 
tational results |14fl2"2"] . We also compare our results to 
those of Kainulainen & Marra |11[ H2] who use a simi- 
lar but slightly different simplified model of large scale 
structure. 



Quantity 


Value 


Mhalo 


1.25 x 1O 12 M 


R s 


0.03 Mpc 


C 


10 



Table I: Parameters of halo with NFW profile 



II. MODEL OF LENSING DUE TO GALAXY 
HALOS AND VOIDS 

A. Galaxy halo profile 

We model the galaxy halos with an NFW profile [ 1 3J 
with a density distribution 



Phaio (r) 



r(r+R s ) 2 



T > CR S 

r ^ CR S 



(2.1) 



Here r is the proper spherical radial coordinate, R s is 
the physical radius which defines the core of the halo 
where most of the mass is concentrated, C the ratio of 
the radius of the halo to the core radius R s , and the 
parameter po is determined by the total mass of the halo. 
The corresponding total halo mass is 



M ha i = AnpoRt log(l + C)- 



C 



l + C 



(2.2) 



For all our simulations we use Mhaio = 1-25 x 1O 12 M , 
R s = 30 kpc and C = 10 [2"lH2"o] . These values determine 
the halo density parameter pq. This completely defines 
our NFW halo model and we list our parameters in Table 
I. In this paper we keep the halo parameters fixed, but it 
would be straightforward to explore other values. 



B. Our void models 

In Swiss cheese models, the Universe contains a net- 
work of spherical, non-overlapping, mass-compensated 
voids. The voids are chosen to be mass compensated 
so that the potential perturbation vanishes outside each 
void. Mass flows outward from the evacuated interior and 
is then trapped on the shell wall. In our previous work, 
[Q , we considered a uniformly underdense interior with a 
(5-function shell on the surface. This model is determined 
by a fixed comoving radius R and by the fraction, /, of 
the total void mass on the shell today. These parameters 
determine the evolution with time of the interior mass 
density and the surface mass density. 

In this paper we generalize the models of pQ to include 
the halo substructure of the voids. We consider two dif- 
ferent idealized models. In the first, each void consists 
of a central, uniformly underdense region surrounded by 
a shell consisting of randomly distributed halos, and in 
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the second, halos are placed randomly both in the inte- 
rior and on the surface. The zero thickness shell is thus 
replaced by halos randomly distributed on the surface of 
the sphere, with the number of halos chosen to match the 
mass of the shell. The number of halos thus evolves with 
time. We call our first model the Swiss Raisin Nougat 
(SRN) model, with "raisins" denoting halos and "nougat" 
the smooth void interior. We call the second model the 
Swiss Raisin Raisin (SRR) model. 

For a given void, we denote by r the physical displace- 
ment from the center of the void at r = 0, and we denote 
by s = ra cx (z) the comoving displacement, where a cx (z) 
is the scale factor of the background ACDM Friedman- 
Robertson- Walker (FRW) cosmology. The quantity that 
determines the lensing magnification is the density per- 
turbation 

(r, z) = p (r, z) - p F RW (z) , (2-3) 

where pfrw i z ) = SHqUm/ (8irGal x (z)) is the back- 
ground FRW density and z is redshift. For the SRN 
model the density perturbation is 

ApsRN (r) = -/ (z) pfwn (z) © (a cx ^void - r) 

N shM (z) 

+ ^2 Phalo (|r - acx^voidfiil) , (2.4) 

i=l 

where the first term is the smoothly distributed interior 
underdensity and the second term is due to halos on the 
surface. Here / (z) is the fraction of the mass of the 
sphere on the surface JT[, Y^oid is the (constant) comoving 
void radius, is the step function, iV s hcii is the number 
of halos on the surface, and is a randomly chosen unit 
vector giving the location of the z-th halo on the surface 
of the sphere. The number of surface halos is 

N shca (z) = f(z)^, (2.5) 

A* halo 

where 

M void = ^7rr v 3 oid a^ x p FRW (2.6) 

is the conserved total void mass. 

For the SRR model, the density perturbation is 

ApsRR (r) = -Pfrw (z) © (aw^void - r) 

Nshell(z) 

+ 2J Phalo (|r - acx^voidnil) 
i=l 

JV coro (z) 

+ Phaio(|r-a ex }wmil)) (2.7) 

z=l 



where the last term represents the halos in the interior. 
Here N COIe (z) = (1 - / (z)) M void /A/ halo is the number 
of interior halos and the vectors rxij are randomly chosen 
in the interior of the unit sphere. 

Now consider a light ray that intersects the void. A 
key role in our computations will be played by the impact 
parameters of the ray with respect to the center of the 
void, and with respect to the centers of the halos. These 
impact parameters will be two dimensional vectors in the 
plane perpendicular to the unperturbed ray. Specifically, 
we introduce a basis of three orthonormal spatial vectors 
ei, e2 and e3 with e3 along the direction of the ray. We 
denote by 

p= V A ^a (2.8) 

A=l, 2 

the comoving impact parameter of the ray with respect 
to the center of the void. We denote by 

b * = H b ? eA = H ° cx [P A ~ r void^f ] e A (2.9) 

A A 

the physical impact parameter of the ray with respect to 
the center of the i-th halo on the surface, where we have 
decomposed the unit vectors as = ^2 A nfeA+nfe^. 
Similar formulae are obtained for the impact parameters 
of the interior halos. 

Even though the SRN and SRR models are highly ide- 
alized, they are more realistic than the void models in our 
previous work pQ . A key feature of our idealized models 
is that they can be evolved in time continuously and very 
simply. Within the context of this highly idealized class 
of models, we study the distribution of magnitude shifts 
relative to what would be found in a smooth cold dark 
matter (CDM) model of the Universe with a cosmological 
constant, A, for different source redshifts. 

It is important to note that our models are not spher- 
ically symmetric, as we break up the shell to form halos. 
We assume that nevertheless the large scale evolution of 
a void is the same as it would be in spherical symmetry. 
We also neglect gravitational clustering of halos on void 
surfaces. Our main aim is to investigate the role of small 
scale clumps in producing magnitude shifts. 

To compute the effects of rays passing through our cos- 
mology, we follow the steps described in Section IIC of 
[lj, with the added halo contributions. Specifically, we 
compute a 4 x 4 matrix for each void, multiply all the 
matrices together, and compute the total magnification 
from the final 4x4 matrix. The explicit expressions for 
the 4x4 matrices in term of line integrals of derivatives 
of the gravitational potential are given in Eqs. (2.18) - 
(2.20) of pQ. We drop all of the integrals over the pro- 
jected Riemann tensor in Eqs. (2.19) of [T] except the 
one in the formula for L A C . We then repeat the com- 
putation N » 1 times to build up the distribution of 
magnifications. 
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Figure 1: Comparison between numerical (points with ranges) 
and analytic results (starred points) for the distribution of 
integrated column depths. 



III. RESULTS FOR A SINGLE HALO 

We now discuss the distribution of magnifications due 
to a single halo. The halo has two distinct regions, the 
core x < R s and the external region R s < x < CR S . Here 
x is the physical distance. To compute the magnification, 
we first compute the lensing convergence, k, analytically. 
For a general density contrast S (s) = Ap (s) / /Ofrw this 
is given by 
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y) 



5{y, z) 



(3.1) 



where y is comoving distance along the ray, ys is the 
comoving distance to the source, a ox (z) = (1 + z)~ , H 
is the Hubble constant, c is the velocity of light, is 
the matter fraction and z = z (y) is redshift. Combining 



the halo profile (2.1) and the second term in Eq. (2.4) 



with Eq. (3.11 gives for the lensing convergence due to 
the halo 
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Here b — |b| is the physical impact parameter 
/ VC 2 -a 2 
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Figure 2: The difference between numerically (k) and ana- 
lytically computed mean convergences (k) as a function of 
number of runs N, for one void, comoving void radius R = 35 
Mpc, and fraction of void mass on the shell today 0.9. Our 
numerical simulations agree with the analytic result to < 1% 
for N = 10 4 . 



tan" 



tan 



C + VC 2 - a + 1 



Here is the step function and a = b 2 /R^. Note that 
a < 1 for 6 < i? s , a > 1 for R s < b < CR S and k = for 
& > Ci? s . 

Finally, the mean of the lensing convergence for a single 
halo is obtained by averaging over the impact parameter 

fCR 3 

/,- = - / n(b)bdb (3.5) 



, i n{b)bdb 
(CR S ) 2 Jo 



8irGa c 



(z) ( „ \ y (ys 
— (PoRs)- 



y) M, 



halo 



ys (Ci? s ) z 

(3.6) 

We define r\ = j p(z)dz to be the integrated column 
density which is proportional to the convergence k. Fig- 
ure 1 is a comparison of log 10 (P(rj)), the logarithm of 
the probability distribution of r\. computed analytically 
(starred points, Eq. (A5l from Appendix A) and the re- 
sults from our code (points associated with the ij bins). 
Within each r/ bin, the mean of the bins agree well with 
analytic results. The width of the bins represents the 
sampling accuracy within those bins, the centers of halos 
are sampled less than the rest of the halos. 

To further assess the accuracy of our numerical results, 
we compute the mean of the distribution for a single halo 
for different numbers of runs (A), and compare this with 
the analytic expression (3.6). We find that the results 



from our numerics agree with the theoretical prediction 
with an accuracy ~ A -1 / 2 as expected. In Figure 2, we 
plot the estimator of the mean as a function of A. We 
see that the accuracy is <1% for A = 10 4 . For the rest 
of the simulations in this paper, we will use A = 10 4 . 
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Quantity 


Value 


Qm 


0.3 




0.7 


Ho 


70 kms^Mpc" 1 


^void 


35 Mpc 


Halo profile 


NFW 


Present fraction of void mass on shell 


0.9 


Fraction of shell mass in halos 


1.0 


Fraction of interior mass in halos 


0.0 



IV. 



Table II: Parameters of SRN model 



RESULTS FOR THE SWISS RAISIN 
NOUGAT MODEL 



The Swiss Raisin Nougat (SRN) model is an idealized 
Swiss Cheese model containing spherical voids with co- 
moving radius F V oid = 35 Mpc. As explained earlier, 
the matter in the interior moves towards the outer edges 
of the void with the evolution of the Universe. For a 
particular void at some redshift, we break up the mass 
on the shell of the void into halos with NFW profiles 
and randomly distribute them on the shell. The mass in 
the interior is smeared smoothly inside the sphere with a 
uniform mass density. The parameters of this model are 
listed in Table II. 

A key change from our void models in [T] is that there 
is no longer a zero thickness shell. One of the issues 
encountered in that model was the logarithmic divergence 
in the variance of the lensing convergence distribution 
due to the zero thickness assumption. Here, however, 
we break up the void surface into halos and the effective 
thickness of the shell is set by the size of these halos which 
acts as a natural cutoff. Hence, the divergence is avoided 
which makes for a more realistic and robust model. 



A. Probability of intersecting a halo 

The expected number of times a light ray hits a halo 
is given by the ratio of the total projected area of all 
the halos in a void to the projected area of the void. The 
expected number of intersections at comoving impact pa- 
rameter p = |p| (comoving distance from the center of the 
void) through the shell at redshift z is 
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Figure 3: The probability distribution of magnitude shifts 
Am for a simulation in the SRN model with Qm = 0.3, with 
one void of radius 35 Mpc at z = 0.45 and sources placed at 
z s = 1, and fraction of void mass on the shell today / = 0.9. 
Top: The probability distribution of magnitude shifts Am 
for a single void for Am positive. Bottom: The probability 
distribution for magnitude shifts Am for a single void for Am 
negative. The total probability for Am < is ~ 0.8. 



Here s is the comoving distance from the center of a void 
and i?haio = CR S is the physical radius of the halo. Av- 
eraging over the impact parameter, p, gives 



JVtat (z) 



f (z) M void (z) R 



halo 



^oid^xW^halo 



(4.3) 



Note that the void radius l^oid is comoving while the 
halo radius i?haio is physical. Both these parameters are 
fixed and do not evolve with time. For example, for a 
void placed at redshift 0, iV s h e ii (0) ~ 0.4. 



B. One void 

In this section we will focus on a single void at redshift 
0.45 and a source at redshift 1. We calculate the ex- 
pected number of halo intersections for a light ray from 
Eq. (4.3). For a 35 Mpc void, using the halo param- 
eters from Section II, we get A s h e ii (0.45) ~ 0.8. For 



N = 10 runs, we keep track of the number of times 
that light rays hit halos and we obtain 8064 instances, 
which agrees well with our prediction. We use the den- 
sity perturbation (2.4) for the SRN model and compute 



the lensing convergence n by summing the result (3.2 1 
over all intersected halos. 
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For the rest of the paper, we will concentrate on the 
distribution of the magnitude shift Am. which is a func- 
tion of lensing convergence k 



2 -5 , 

Am = m 

In 10 



(i-k) -r 



(4.4) 



Here 7 is the shear, which we will discuss in the next sec- 
tion. Figure 3 shows the probability distribution P (Am) 
of magnitude shifts Am we obtain for a single void placed 
at z = 0.45 (without shear). 

A notable feature of this distribution is that it is bi- 
modal, with peaks at both positive and negative Am (We 
plot separately the distribution for positive Am and for 
negative Am, since the relevant scales for these two re- 
gions of the probability distribution are very different). 
The peak at positive Am is predominantly due to rays 
that do not intersect any halos, and are demagnified by 
their passage through the underdense void interior. The 
peak at negative Am is predominantly due to rays which 
intersect one or more halos and are consequently magni- 
fied. 

We can also compute the mean of each of these distri- 
butions and compare them with analytical expressions. 
The means for Am > 0, (Am) + and for Am < 0, (Am)_ 
are defined to be 



(At 



AmP (Am) d (Am) , (4.5) 



(Am) = / AmP (Am) d (Am) . 

J — OO 



(4.6) 



We decompose the full distribution of magnification as a 
sum 



P (Am) = J2 V nPn (Am) , 



(4.7) 



where P„ is the probability of n halo intersections and 
P n (Am) d (Am) is the probability of having a magnitude 
shift between Am and Am + d (Am) given that there are 
n intersections. 

The analytic expression for the magnitude shift for zero 
halo intersections is |T] 



Am = 2.5 log 



10 



(1 - K h 



r)' 



(4.8) 



where 



c ysa e x(z) v 



Y 2 

void 



(4.9) 

Here / (z) is the fraction of the mass of the void on the 
shell, y is the comoving distance to the void, ys is co- 
moving distance to the source, a ex (z) is the scale factor 
and p is the comoving impact parameter. Note that for 
this one void case, essentially all the negative Am con- 
tributions are due to intersections with one halo and the 



positive Am contributions are due to the void interior 
(i.e., no halo intersections). In this approximation, Eqs. 
(4.5) and (4.6 1 reduce to 



and 



/OO 

AmP (Am) d (Am) (4.10) 
-OO 

/OO 

AmPi (Am) d (Am) . (4.11) 
-OO 



We can compute Vq and V\ from Eq. (4.3) assuming 
P0+P1 = 1 for the one void case, obtaining V\ — 0.8 and 
Vq = 0.2 which matches with our simulations. The nu- 
merically computed means (4.5) and (4.6) of the magni- 
fied and demagnified distributions agrees with their cor- 
responding halo and void interior theoretical values [com- 
puted from Eqs. Q - Q, fO} & to - 0.5%. 
We also numerica 
gence obtaining 
viation 2 x 10 



Uy compute the mean lensing conver- 
-5 x 10 -4 magnitudes with standard de- 
" 3 magnitudes. Thus the mean is consistent 
with zero as we would expect from a general theorem. 



C. Shear 

So far in our analysis we have not included shear. We 
can include it as follows. The matrix L A C for the j-th 
void defined in Eq. (2.15c) of [1] is 



L AB = -2 dy 



(4.12) 



where 5$ is the potential perturbation, y is comoving 
distance, the derivatives are with respect to comoving 
coordinates, and the integral is taken over just the j-th 
void. We decompose Lac into a trace part and a trace 
free part to obtain 



Lab = — [ft>j Sab + IjAb] , 



(4.13) 



where Wj — yj (ys — Vj) /ysi Vj is the comoving distance 
to the j-th void, Kj is the lensing convergence we com- 
puted previously [Eqs. (3.2 1 - (3.4 1 and (4.9 1], and the 



matrix 7 is traceless. We compute the potential per- 



turbation from the density perturbations (2.3) and (2.4 1, 
and insert into Eqs. (4.12 I and (4.13 I to obtain the shear 



term jjAB which will be of the form jjeAB on a suitable 
choice of basis. For a single NFW halo the potential per- 
turbation is 



8<f> (r) = -AnGpaRl 



/Js In I 1 , 



1 



1 + C 

(4.14) 

The shear due to the void interior and the halos on the 
shell is 
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Here p = pp A ep^ is again the comoving impact parameter 
to the void, P — p 2 /Y 2 oid and = bibfe^ is the physical 
impact parameter of the light ray to the z-th halo. The 
first term in Eq. (4.17) is the point mass contribution 
of the halos. The second and third terms in Eq. (4.17) 
are non zero only for intersected halos. Evaluating the 
integrals using the NFW profile $2.1 s. \w Mud tluii i in- 
intersected halo contribution 
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Here for a < 1 
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1 (2 - a) VC 2 - a , 
lh= a 2( a -l )( C + l) +ln 



(3a - 2) 
\l-af 2 



tanh 1 



c + Vc 2 





' ll-y/a' 


^tanh 1 







C + l + ^/C 2 ~ 



a 



(4.19) 



and for 1 < a < C 2 



1 (2 - a) VC 2 - a , 
lh= a 2(a-l)(C + l) +ln 



^- 2 ) ^tan- 



(a-l) 



3/2 




tan 



C + l + \/C2 



(4.20) 



where a = bf/R 



From Eqs. ( |4.18| - ( |4.20| , 7 2 - 1/6 4 . Therefore con- 
tributions from shear are heavily suppressed. In our nu- 
merical analysis, the standard deviation changes by less 
than 3% if shear is neglected. 



Figure 4: The probability distributions of magnitude shifts 
Am for the SRN model with sources at redshifts of z s =0.5 
(dashed), z s — 1.0 (dot-dashed) and z s — 1.5 (solid), for 
comoving voids of radius R = 35 Mpc with 90% of the void 
mass on the shell today. The horizontal lines are the 25% 
(top), 50% (middle) and 75% (bottom) quartiles about the 
peak of the distribution. 



D. Qualitative features of magnification 
distributions 



With the accuracy of our method tested, we now ex- 
plore the magnification distributions in more general sit- 
uations with many voids, distributed along the line of 
sight with random impact parameters according to the 
algorithm discussed in [Q. For example, for sources at 
rcdshift z s = 1, there are 47 voids of comoving radius 
^void = 35 Mpc along the line of sight. We follow steps 
1 to 8 of Section IIC of our previous paper [T] , but with 
the modification that the matrices J, K, L and M now 
incorporate the effects of the halo substructure of the 
shell. 

In Figure 4, we plot the log of the magnification dis- 
tribution for z s — 0.5, 1.0 and 1.5. In our SRN model, 
we have voids with randomly distributed halos on their 
surface and a smooth interior. We denote by V n the 
probability of having n halo intersections. The whole 
probability distribution can be decomposed into a sum 
of probability distribution for different numbers of halo 
intersections, like in Eq. ([47}. We plot the 25%, 50% 



and 75% quartiles of the distributions as horizontal lines 
(top, middle and bottom respectively). For high red- 
shifts, most of the probability is concentrated in the de- 
magnified areas where the rays hit only a few halos or 
simply pass through without hitting any. 

Note that the total probability in the tail on the mag- 
nification side (Am < 0) increases with redshift, because 
of the increased probability of hitting halos at higher red- 
shifts. For example, at Am — —0.2, we would expect the 
probability density for z s — 1.5 to be roughly 2-3 times 
as large as the probability density for z s — 0.5 because 
the number of voids that rays have to pass through in the 
former case is 62 where as for the latter it is 27. In ad- 
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dition, rays at high redshifts have more close encounters 
with halos that generate shear. 

Figure 5 shows the standard deviation of the distri- 
bution, function of redshift of the source, z s . 
This standard deviation for voids and halos is ~ 3 times 
larger at z s = 1 than that for a model with mass com- 
pensated voids with a shell thickness of 1 Mpc and no 
halos [1]. We note that most of the contribution to the 
standard deviation come from rays that intersect halos. 
Also, the standard deviation we compute agrees well with 
that computed using other methods. For example, our 
standard deviation for z s = 1.5 is cr m = 0.072, which 
agrees to within 20% with the standard deviation of the 
distribution shown in Figure 1 of Ref. [5]. We compare 
our results to those obtained using another method in- 
troduced in Refs. IT2~] in the next subsection. 

In Appendix B we derive the following approximate 
result for the standard deviation: 



log (10) 



3 



(HoY^f^H^w 



n 2 



E 

3 



(-ffo^void/j 



) H w j 



(4.21) 



Here g h is a dimensionless parameter which represents 
the contribution from halos whose detailed form is given 
by Eq. (B21), Y vo id is the comoving radius of the voids, 
fj = f (zj) is the fraction of the mass of j-th void on 
its surface and Oj = a (zj) is the scale factor. The result 
(4.21) assumes statistical independence of halos within 
voids and also of voids from one another and neglects 
lensing shear. There are two main qualitative features of 
the result (4.21). First, the contribution to the standard 
deviation due to the halos [the first term in Eq. (4.21)] 
depends primarily on their gravitational potential, and 
the contribution due to voids [the second term in Eq. 
(4.21)] depends primarily on the size of the underdense 
core. Second, the halo contribution is bigger than the 
interior contribution and hence the standard deviation 
is dominated by halos. For example, using the above 
expression, the ratio of the contribution due to the halos 
to the contribution due to the core is ~ 100 for z s = 1 
and for the void and halo parameters defined in Table II. 

We discuss further the analytic calculation of standard 
deviation without shear in Appendix B. Our numerical 
results agree with these approximate analytic predictions 
to within ~ 20%. 



E. Redshift dependence of mean and mode of 
magnitude shift 

While the mean of lensing convergence vanishes, the 
mean magnitude shift does not, because magnitude shift 
is a nonlinear function of k, defined in Eq. (4.4 1. Figure 
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Figure 5: Redshift dependence of standard deviation of dis- 
tribution of magnitude shifts, for comoving voids of radius R 
= 35 Mpc with 90% of the void mass on the shell today. The 
crosses are analytic results. 



6a shows the mean /x m of the distribution of magnifica- 
tion shifts Am, which increases with redshift as oc <j 2 n . 
This is the expected theoretical behavior: for small values 
of k and ignoring shear, we can approximate Eq 
as 



(4.41 



In 10 



In |(1 - k) 



5 

In 10 



1 , 

— K 

2 



(4.21) 



The mean magnitude shift is then proportional to the 
mean of the square of k as the mean of k is vanishing, 



Mr. 



2.5 
In 10 



(4.22) 



The standard deviation, from Eq. (3.4 1 in [T] simplifies 
to 



In 10 



V¥ 2 



(4.23) 



and so fi m ~ — 0.23cr T 2 n which agrees with our numerical 
results to within ~ 10%. 

On average there is a small overall magnification of 
light beams. Figure 6b shows the mode i/ m , the location 
of the maximum of the PDF, which also increases with 
redshift. The modes of the magnification to redshift 1.5 
are positive because an overall demagnification occurs for 
most of the light rays as they pass through the interior 
of the voids while hitting halos. Note that the modes are 
larger than the corresponding means. 

In realistic surveys, one can expect to find only a few 
standard candle sources for every redshift or every red- 
shift bin. This severely constrains the accuracy of the 
cosmological parameters we can infer from such observa- 
tions. To illustrate the extent of lensing degradation in 
measuring cosmological parameters, we pick 200 sets of 
randomly placed 10 or 100 sources at z s — 1.5. We find 
the mean of each of these samples in the set and plot the 
resulting distribution of means in Figure 7 (top). The 
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Figure 6: Plot of mean (top) and mode (bottom) of the dis- 
tribution as a function of redshift in SRN. The mode takes 
on increasingly positive values with redshift and the mean is 
increasingly negative with redshift. 



mean of the means for the 10 sources case is 0.011 mag- 
nitudes and the standard deviation of the means is 0.028 
magnitudes. The respective numbers for the 100 sources 
case are 0.0004 and 0.011 magnitudes. A change in cos- 
mological parameters by 1% implies a change in Am of 
0.015 magnitudes. Thus for data acquired from surveys, 
the lensing degradation is quite a significant effect, al- 
though it can be mitigated by increasing the number of 
sources. This is also seen in Figure 7 (bottom) where we 
plot the range of standard deviation for samples of differ- 
ent sizes. For a large enough sample, the bias in magnifi- 
cation can be accurately taken into account. This effect 
is studied in [5] which shows that lensing degradation ef- 
fectively decreases the number of useful supernovae by a 
factor of 3 at source redshift 1.5. 

Our work is broadly consistent with other work, |lll 
H2l fT4"H2"2"] in this area. A similar computational method 
has been developed by Kainulaincn & Marra, Refs. [11, 
12 . Their model consists of filaments and halos of various 
sizes, where the mass fraction in filaments is 0.5 and the 
rest is distributed in halos. To compare with their results, 
we use the SRN model and choose parameters to match 
their cosmology, i.e., 17m = 0.25, z s — 1.5, H Q = 73 
kms _1 Mpc _1 and / = 0.5. We do not include shear 
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Figure 7: Top: Plot of distribution of means of 200 samples of 
10 (dashed) and 100 (solid) sources each at redshift z s — 1.5. 
The mean of the distribution of means is 0.011 magnitudes 
and the standard deviation is 0.028 magnitudes for the sample 
of 10. The respective numbers for the sample of 100 sources 
is 0.0004 and 0.011, showing that more sources reduces de- 
magnification bias. Bottom: Plot of the range of standard 
deviation for samples of different sizes showing convergence 
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for this comparison as it is neglected in their analysis. 
Our magnification PDF qualitatively agrees with that of 
Kainulainen & Marra as shown in Figure 8. 



V. RESULTS FOR SWISS RAISIN RAISIN 
MODEL 

In our Swiss Raisin Raisin (SRR) model, in addition to 
replacing the smooth surface density on the shell with a 
collection of halos, the mass in the interior is also broken 
up into NFW halos with the same parameters as before. 
The parameters of the SRR model are listed in Table III. 



A. A single void 

In this section we focus on a single void with no shear. 
Again, we use Eqs. dO), pH} - pll) and (|4.8h to com- 



pute the magnifications. For a 35 Mpc void and using 
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Figure 8: The probability distributions of magnitude shifts 
Am for an SRN model (dashed lines) with sources at redshift 
z s = 1.5, comoving voids of radius 35 Mpc, with 50% of the 
void mass in halos today, = 0.25, Q,a = 0.75, Ho = 73 
kms _1 Mpc~ 1 and with no shear, compared to the results in 
Figure 5 (we reproduced the plot by picking points from their 
figure) of the model in Kainulainen and Marra [12, where they 
have 50% of mass in halos and all other parameters same as 
ours. The two distributions are qualitatively similar. 



Quantity 


Value 




0.3 


n A 


0.7 


H 


70 kms^Mpc" 1 


^void 


35 Mpc 


Halo profile 


NFW 


Present fraction of void mass on shell 


0.9 


Fraction of shell mass in halos 


1.0 


Fraction of interior mass in halos 


1.0 



Table III: Parameters of SRR model 



the same halo parameters as in Section II, the intersec- 
tion probability remains about the same as in the pre- 
vious model. The change due to the addition of a few 
halos in the vast interior region is negligible. The ex- 
pected number of halo intersections Ni nt (z) is given by 
Eq. (4.3) with the shell mass replaced by the total void 



mass A'L 



void 



N int (Z) 



M v 



mR 2 



halo 



r void a cx ( Z ) M halo 



(5.1) 



Using the same parameters as in the one void case in the 
SRN model we obtain N int (0.45) ~ 0.88. We find that 
light rays hit halos 8720 times for N = 10 4 runs which 
translates to an intersection probability of ~ 88% and 
agrees well with Eq. (5.1 1. 



Again we find a bimodal distribution of magnitude 
shifts. In Figure 9, this bimodal distribution is superim- 
posed on the SRN plots from Figure 3. Since the density 
contrast in the interior is increased by ~ 10% (if halos 
are not hit) compared to our SRN model, the demagni- 
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Figure 9: The probability distribution of magnitude shifts 
Am for a simulation in the SRR model (dashed line) super- 
imposed on the corresponding probability distribution in the 
SRN model (solid line) with Q m = 0.3, with one void of radius 
35 Mpc and fraction of void mass on the shell today / = 0.9. 
Note that the demagnified part, (top), is shifted because there 
is an increase in the density contrast in the interior of voids, 
while the magnified part, (bottom), is mostly unchanged. 



fied distribution shifts towards the right. Figure 9a shows 
this shift in underdense part of the distribution. Figure 
9b is the Am < part of the distribution. This accounts 
for roughly 88% of the total distribution and it is similar 
to the distribution in Figure 3b. 



B. Redshift dependence of distributions 

Next we explore the bias (i.e., the mean of the dis- 
tribution) due to halos and voids for sources at various 
redshifts, namely, for z s — 0.5, 1.0 and 1.5. To compute 
the shear in this case, for the void contribution we use Eq. 

1, and in Eqs. Ell} 



(4.201 



( gg) ) but with / (z) 
sum over both the surface and interior halos. The mag- 
nifications shown in Figure 10 are predominantly due to 
halos while the mostly empty interior has a demagnify- 
ing effect. Due to the increased underdensity inside the 
void, the magnitudes shift towards demagnification. The 
standard deviation for z s — 1.5 is a m — 0.072. Again we 
note that the modes are larger than the means and also 
shift towards the demagnification end of the plot with 
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Figure 10: The probability distributions of magnitude shifts 
Am for the SRR model with sources at redshifts of z = 0.5 
(dashed), z = 1.0 (dot dashed) and z = 1.5 (solid), for voids 
of comoving radius R = 35 Mpc, with 90% of the void mass 
on the shell today. We see the same qualitative features as in 
SRN for the corresponding redshift but all the distributions 
are shifted towards demagnification. The horizontal lines are 
the 25% (top), 50% (middle) and 75% (bottom) quartiles of 
the distribution from the peaks. 



increasing redshift. The tails of these distributions are 
similar to the ones obtained in the SRN model in Figure 
4. 

Next we consider the mean magnitude shift, f^m^ 

and 

its mode, v m , in the two models. The key feature here 
is that the underdense interior is more prominent in the 
SRR model. We expect that the mean magnitude shift 
and its mode should shift and the difference in the means 
should increase with redshift. In Figure 11, we plot the 
means and modes for the two models and we show that 
Hm for SRR is - 10 - 20% greater than that for SRN at 
z s = 1. 

These two models are interesting because they repre- 
sent the two possible extremes of the matter distribution 
in voids - one where the matter is smoothly distributed 
with no structure and another with only chunky NFW 
halos. In reality, the underdense region will be composed 
of both halos and an ambient intergalactic medium. By 
studying the completely smooth interior case (SRN) and 
the completely granular interior case (SRR), we expect 
to bracket the true distribution. 



C. Effect of large scale structure 

Previous studies in the literature have modeled the 
magnification effects of only voids [TJ UJ [TU] or only halos 
[iiOOQT], while other studies have considered very specific 
models with a particular distribution of filaments and ha- 
los [12]. In our work, we present two models that incor- 
porate the effects of both halos and voids. We have con- 
sidered cosmological models where at large scales matter 
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Figure 11: The mean and mode of magnification shift for the 
two models - Points: SRN; Crossed points: SRR. Top: The 
mean of magnification shift, fi m , for the two models. We see 
that the difference in the means increases with redshift and at 
z = 1 and it is ~ 10%. Bottom: The mode of magnification 
shift, u m , for the two models. We see that the difference in the 
means increases with redshift and at z = 1 and it is ~ 20%. 



evolves to cluster on the edges of spherical voids. 

One limit of our model is the case where there are no 
halos on the surface and the interior is composed entirely 
of halos. This corresponds to the / = limit of the SRR 
model, a "no void limit". In Figure 12, we compare the 
magnification distribution we obtain for this configura- 
tion for z s = 1.5 to the corresponding distribution in the 
SRR and SRN models, for comoving voids of radius R = 
35 Mpc with 90% of the void mass on the shell today. 
From Eq. (5.1 1, the expected number of halo intersec- 



tions in the SRR model is independent of the radial dis- 
tribution of the halos and hence the SRR and no void 
distributions should be similar, as observed. However, in 
the SRN model the number of halo intersections is lower 
due to the smaller number of halos. Hence the distribu- 
tion is shifted towards negative magnification shifts. 
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Figure 12: The probability distributions of magnitude shifts 
Am for simulations in SRR (solid), the SRN model (dot 
dashed), for comoving voids of radius R = 35 Mpc with 
/ = 0.9 today. Also shown is the model with uniformly dis- 
tributed halos everywhere and with no voids (dashed - this 
is the / = limit of the SRR model) with source at red- 
shifts of z = 1.5. The distributions are similar in the SRR 
models because the expected number of halos intersected is 
independent of of the radial distribution of halos, where as 
the distribution in the SRN model is qualitatively similar but 
is shifted towards negative magnifications. 

VI. CONCLUSIONS 

In this paper, we presented two simple models to study 
the effects of both voids and halos on distance modulus 
shifts due to gravitational lensing. Our results may be 
useful for future surveys that gather data on luminos- 
ity distances to various astronomical sources to constrain 
properties of the source of cosmic acceleration. The core 
of our model is constructed by considering a ACDM Swiss 
cheese cosmology with mass compensated, randomly lo- 
cated voids, while our small scale halo structures are non- 
evolving and chosen to be all the same size and with an 
NFW matter profile. 

We used an algorithm, described in |T|, to compute the 
probability distributions of distance modulus shifts. The 
mean dispersion of the magnitude shift due to gravita- 
tional lensing due to voids and halos is ~ 3 times larger 
than due to voids alone with a shell thickness, [T]; the dis- 
persion a m due to 35 Mpc voids and halos for sources at 
z s = 1.5 is a m = 0.065 — 0.072 (depending on the model). 
The mean magnitude shift due to voids and halos is of or- 
der 5m = —0.0010 to —0.0013 (depending on the model). 
These values of a m imply that large scale structure must 
be accounted for in using luminosity distance determi- 
nations for estimating precise values of cosmological pa- 
rameters, such as those characterizing the dark energy 
equation of state. 

We studied the distribution of magnitude shifts for 
three different source redshifts in each of our models. 
The qualitative dependence on redshift is similar to that 
of the previous void-only models [1 . We find that the 



voids do not significantly change the variance but do sig- 
nificantly change the demagnification tail and the mode. 
The mode lies on the demagnification side and the vari- 
ance is largely due to halo intersections. The scale of 
voids is unimportant and the only discernable effect in 
the mode is seen when the void interior is smoothly dis- 
tributed matter. As a result, our models bracket the 
range of possibilities of magnifications. 

Our simple and easily tunable model for void and 
galaxy halo lensing can be used as a starting point to 
study more complicated effects. For example, one can 
use various algorithms to generate realizations of distri- 
butions of non-overlapping spheres in three dimensional 
space. Given such a realization one could use the algo- 
rithm of this paper to study correlations between magni- 
fications along rays with small angular separations, which 
would be relevant to future small beam surveys [23]. Fi- 
nally, our model is in general agreement with other sim- 
plified lensing models in the literature that focus on lens- 
ing due to both halos and larger scale structures, [121 HU- 

Our results for <j m in the SRN void model are repre- 
sented within about 20% by an analytic model presented 
in detail in Appendix B. This model ascribes the magni- 
tude shift entirely to the fluctuations in light beam con- 
vergence that results from passage through underdense 
cores and overdense halos; thus it ignores the contribu- 
tion from shear, which we have found to be relatively 
small empirically. The final result for is a sum of 
these two contributions. [See Eqs. (B20) and (B22)]. 
Although our simulations assumed a single halo mass 
M/j, radius and concentration C, and a single void 
radius, the analytic model allows distributions for these 
key quantities. The contribution from halos is propor- 
tional to a suitably weighted mean of Mh^2{Ch) / 
where ^(Cft) is defined in Eq. (B16) and is displayed 
in Fig. 13. The contribution from the underdense void 
cores is proportional to the mean void radius. Typically, 
the contribution to from halos is much larger than 
the contribution from void cores so ct^ is larger for more 
massive or more compact halos. 

We have seen that the results of our simulations depend 
on whether the underdense core consists of smoothly dis- 
tributed dark matter (SRN) or is itself clumped into halos 
(SRR). In the extreme case in which the underdense core 
is entirely made of halos, the results do not depend on 
the core density, and is equivalent to the SR model that 
consists of halos distributed randomly within a void. 

Although the analytic model was only developed for 
the SRN model, it could also be applied to the SRR 
model with any prescription for the fraction of the mass 
of the underdense core that is clumped into halos. For 
example, in the extreme case of total clumping we can 
use the analytic models in Appendix B with the substi- 
tution / = 1 for all z in all expressions derived there. 
For intermediate prescription for the fraction of 

the underdense core that remains smooth rather than 
clumped would be needed. 
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APPENDIX A: COLUMN DEPTH 

In this appendix, we describe how we calculate the 
column density encountered by rays passing through ha- 
los in the SRN model. In our model for the voids, we 
break up the bounding shell of mass into halos with mass 
-Whaio- As a light ray passes through one of these ha- 
los, the beam will acquire some integrated column depth, 
r\ = J p(z) dz. where the random variable rj depends on 
the impact parameter of the beam with respect to the 
halo center, p (z) is the density profile of the halo and 
y is the physical coordinate. The maximum value is for 
a beam going right through the center, and diverges for 
our NFW profile pH] 



We use an NFW profile 2.1 for the matter distribution 
in halos, and the column depth is 



rv /(Cfl s ) 2 -b 2 dz 

V (b) = (ZpoRl) / - - % - 2 

Jo v'z 2 + b 2 (y/z 2 + b 2 + R s ) 

(Al) 

where we have changed from radial (in Eq. ( |2.1| ) to 
Cartesian coordinates. Here the physical impact param- 
eter is b and C is the ratio of the radius of the halo to its 
core radius. This reduces to the sum of two contributions 

?7haio (b) = 77cor c e (R s - b) + Vout e (b - R s ) 9 {CR S - b) 

(A2) 

where the relationship between the column depth rj (b) 
and the corresponding lcnsing convergence, k (6), is, from 
Eq. (3.1), 



rj(b) 



( 4irGa ex (z) y (y$ - y) 
V c 2 y s 



«(6). 



(A3) 



The lensing convergences is listed in Eqs. ( 3.3 ) and ( 3.4 1 



Again, the contribution from outside a radius of CR S is 
zero. The mean column density of halos is defined as 



?7halo 



alo 



w(CR s 



(A4) 



For the halos, even more important that rj is drj/da, 
where a = b 2 /R 2 ; this is because the probability dis- 
tribution for a is da/C 2 , and therefore the probability 
distribution for rj for a single halo is 



-Phalo (r?) 



1 



C 2 \dri halo /da\ ' 
This is the quantity plotted in Figure 1. 



(A5) 



APPENDIX B: ANALYTIC ESTIMATE OF 
STANDARD DEVIATION 

In this appendix, we derive analytic results for the 
standard deviation of magnifications. Consider the mean 
of the contribution to the convergence (3.1) from the un- 
derdense core of the j-th void, which we will denote by 
K c j. We find 



3H Qm 



fjVjjys - v 3 ) 

ajVs 



x 2 



Y 2 - p 2 

3 1 3 



(Bl) 



where Yj is the comoving radius of the void, aj is the 
scale factor and fj is the fraction of the total void mass 
on the surface at redshift Zj . After averaging over impact 
paramters we obtain 



-2n M H Y ^ Hoyjiys - Vj \ (B2) 
a ys 



and therefore the net expected convergence from voids is 



fjYjHoyjjys ~ Vj) 
ajys 



(B3) 

Assuming a typical radius Yj ~ ^void, there are about 
(-ffo^void) _1 terms in the sum, and consequently the over- 
all average is ~ &m- I n the limit that there are many 
voids along a given line of sight, we can replace the sum 
by an integral. The number of voids per interval dy of 
comoving distance is dy/2R, and therefore if we define 
£ = H y we find 



3 
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(B4) 



where £5 = Hoys- Equation (B4) does not depend on 
any void properties (apart from the value of / today) 
and remains valid if there is a distribution of void sizes, 
for example. 

On average, the contribution k/j from halos to the lens- 
ing convergence must cancel the contribution (B4) from 
voids, i.e., (« h ) = - (k c ). 

We assume statistical independence of halos within 
voids from the core (i.e., true if halo radii are small) 
and also of voids from one another. The overall vari- 
ance is a sum of individual halo and core variances. As 
in our simulations, we neglect clustering of halos, which 
would introduce correlations among them, and assume 
that dark matter is confined to the halos and underdense 
core. These assumptions could be relaxed in a more so- 
phisticated model. We see that 
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since the averages for i =/= j vanish. The sum is < 
(HqY)- 1 and therefore < fl 2 M H Y. 

Let us now consider halos residing in void j. We in- 
clude the possibility that there are different types of halo 
with different properties, and label the types by a. For a 
given halo i of type a passed through by the line of sight 
at physical impact parameter b a ^ relative to its center, 
the contribution to Kh is 



8-KGajyjiys - yj) 



Us 



dr r p a (r) 



r 2 - b 2 - 

Oi,l 
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where Rh, a is the physical halo radius, and p a (r) is the 
physical density within the halo. The average over im- 
pact parameters b a i is 



R l 



db a {b a: 



dr r p a (r) 



r 2 - b 2 ■ 



the distribution of halos within the void but depends on 
the assumption that all halos are associated with voids. 
Let us suppose that 



4GMh, a a a yj(y s - yj) I b 2 a ,i C h, a 
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where Ch, a is dimensionless. This form assumes that the 
density profile for halo type a has one scale parameter, 
Rh.a/Ch.a, although it does not necessarily assume that 
the density profiles are the same for all halos either in 
form or in the parameter Ch a - We do assume that 



dxF a (xCl a ) = l = 



1 



CI 



dqF a (q), (B13) 



independent of the value of Ch.a- With this normaliza- 
tion, if we let p a (r) = po, a Pa(rC h , a / Rh,a) we find 
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(B7) 



where Mh, a is the total halo mass. Therefore the average 
over impact parameters through a given halo is 



AGMh^ajyjjys - yj) 
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If all of the halos reside in the voids (i.e., none in the FRW 
exterior), then the expected number of intersections of 
the light path with a halo is iV QJ i?| a /a 2 Y 2 , where N a j 
is the expected number of halos of "type a" in the void, 
and we take account of the fact that Rh, a is a physical 
radius. Then we get a total contribution from halos per 
mass-compensated void equal to 

\- N <*,j R l, a ( K (c*,i),j) 4Gy 3 (ys - Vj) \- 

2^ 7^72 - T172 N <*,j M h,a 

a j 3 ' ' : 



(B9) 

the sum is the total mass in halos, which must compen- 
sate the underdensity, so 
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and therefore 

^aj^aka) _ 2f j H${l M Y j y j (y s - Vj ) 
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Therefore the average per mass-compensated void cancels 
as expected. This cancellation is actually independent of 
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Equation (B12) with the normalization (B13) leads to 



Eq. (Bill as expected, but it also implies contributions 



from each halo to the variance given by 



2 / \ 2 

K (a,i),j \ K (a,i),j) 
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Averaging over impact parameters implies 



( K (a.i),j) ~ ( K (a,i),j) 
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where the function ^2(Ch, a ) for NFW profiles is plotted 
in Figure 13. 

We can use these results to determine the expected 
contribution of halos in a given mass-compensated void 
to the variance. Since total mass in halos in void j is 



15 



15 - 
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the dispersion. There is also a hefty numerical factor 
8^2(Ch. a ) ~ 50 favoring the contribution from halos. 

If Va,i — Va is actually independent of i then we can 
factor out the sum over a in the first term of Eq. ( B20 1 : 
define a dimensionless parameter 
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0.175(M h , a % 2 (C h , a )/Ri a ) 
h x 10 12 M Q /(300kpc) 2 x 6.84. 



(B21) 
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Figure 13: Plot of the behavior of \&2(Ch, Q ) as a function of 
the concentration parameter Ch, a of NFW halos. 



We rewrite Eq. flB20| as 



M- \ ^ 
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fj^MH^R 3 - J2G , if the fraction of this total in halos of 
type a is rj a ^, then expected number of type a is 



N ■ = 



2GM Ka 

The expected number of intersections with type a is 
ri a! ifjClMH$YjRl 



(B17) 



2a]GM hi , 
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The contribution from halos in mass-compensated void j 
to the total variance is therefore 

8/j GCImHqYj y 2 (y s - y 3 ) 2 r, ati M h ^ 2 {C h , a ) 



E 



(B19) 

Summing over all mass-compensated voids out to the 
source we get 
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n h,ot 
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Equation ( |B20 1 shows that the overall variance depends 
on halo properties primarily via the typical value of 
GMf l a /R 2 la , the gravitational acceleration characteris- 
tic of the outer regions of the halo. When this is large 
compared with the mean gravitational acceleration of 
the void as a whole, fjflMHoRi/a 2 , the halos dominate 
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The first term in Eq.( B22| ) is ~ g h and the second is 
- fl M Ho(Yj) « 3.5 x 10- 3 /i(Yj)/35Mpc. Eq. p21 | 
suggests that halos dominate. 

Now all of the parameters in Eq. ( B22 1 can be varied 
over distributions. To keep things as simple as possi- 
ble, let us assume that fj is simply a function of redshift; 
that is, neglect the possible dependence of / on the size of 
mass compensated voids. As a further simplification, we 
turn the sums into integrals. If we had a single void ra- 
dius R the sums would turn into integrals by noting that 
there are dy/2R voids per range dy: for a distribution of 
void sizes we can use this substitution with Ri — > (R) 
in the second sum. If we also define £ = Hoy then with 
these simplifications 
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The relationship between the variance of k and the vari- 
ance of magnifications a m , for small deviations, is ap- 
proximated by 



(B24) 



log (10) 
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